Metastable states of a gas of dipolar bosons in a 2D optical lattice 
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We investigate the physics of dipolar bosons in a two dimensional optical lattice. It is known 
that due to the long-range character of dipole-dipole interaction, the ground state phase diagram 
of a gas of dipolar bosons in an optical lattice presents novel quantum phases, like checkerboard 
and supersolid phases. In this paper, we consider the properties of the system beyond its ground 
state, finding that it is characterised by a multitude of almost degenerate metastable states, often 
competing with the ground state. This makes dipolar bosons in a lattice similar to a disordered 
system and opens possibilities of using them as quantum memories. 



The effect of long range interaction on the quantum 
phases of ultracold gases in optical lattices has been re- 
cently investigated in the literature [3, 0]. Theoretical 
studies have pointed out that novel quantum phases, like 
supersolid and checkerboard phases, arise as soon as the 
interaction potential involves at least one nearest neigh- 
bour. This issue has recently become of primary impor- 
tance, because, since the achievement of Bose-Einstein 
condensation of dipolar Chromium atoms [H and with 
the progress in cooling and trapping of dipolar molecules 
[H, these systems are becoming to be at experimental 
reach. 

As mentioned above, the phase diagram for non-zero 
range interactions presents two main kinds of phases: (i) 
superfluid supersolid phases, where large atom-number 
fluctuations are found at each lattice site and an order 
parameter different from zero characterises the system 
M [fit; (ii) insulating checkerboard phases, where num- 
ber fluctuations are absent and a well defined number of 
atoms is found in each lattice site. These two phases dif- 
fer from the usual superfluid and Mott insulating phases, 
because they present modulated patterns in the density 
and in the order parameter (when different from zero). 
Usually these patterns consist in regular distributions of 
atoms in the lattice sites, and are characterised by a fill- 
ing factor (average number of atoms per site) which is in 
general not integer, even in the insulating phases. 

The clear observation of such phases is a very impor- 
tant experimental challenge. Although the observation of 
a possible supersolid phase in 4 He has been reported Q , 
another proof of the existence of such quantum phases is 
desired. We think that it could be eventually obtained 
with samples of ultracold atomic gases in optical lattices 
in the presence of long range interaction [8( . 

In this work, we focus our attention on the insulat- 
ing states which are found in the low-tunneling region 
of the phase diagram. We point out the existence of 
metastable states in the system. Contrary to the stan- 
dard on-site Bose-Hubbard model, we find that beyond 
the ground state, there exist a huge amount of config- 
urations ("classical" distribution of atoms in the lattice 
sites), which have higher energy but result stable against 
tunneling. For small lattices (4 x 4) and periodic bound- 



ary conditions, we analyse all possible existing configu- 
rations and we find that for given chemical potential and 
tunneling parameter there can be as much as hundreds of 
metastable configurations. We then generalize our results 
to large lattice sizes (typically up to 20 x 20) , comparable 
to the ones found in realistic experimental situations. 

Characterising systematically the metastable states in 
terms of their stability against perturbations and their 
capability of being approached in the time evolution of 
the system in the presence of dissipation is important 
in view of the possible application of those systems as 
quantum memories. To our knowledge, the existence of 
metastable states has not been discussed in the litera- 
ture on extended Bose-Hubbard models so far Our 
results are based on a mean-field description of the sys- 
tem, but we believe that the existence of the metastable 
states should be confirmed by treatments beyond mean- 
field which are sensitive to local minima of the energy. 

We consider a two-dimensional (2D) gas of dipolar 
bosons in the presence of a 2D optical lattice, and an ex- 
tra confinement in the perpendicular direction [111 ]. We 
assume a single component gas of bosons (i.e. spin, or 
pseudo-spin, polarised) [l2l ]. Our system is well described 
by the extended Bose-Hubbard Hamiltonian 
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where J is the tunneling parameter, U the on-site interac- 
tion, U? the components of the dipole-dipole interaction 
at different relative distances, and \i the chemical poten- 
tial which fixes the average atomic density. The notation 
(ij) represents nearest neighbours and ((ij))? represents 
neighbours at distance £. 

The on-site interaction is given by two contri- 
butions: one is arising from the s-wave scattering 
U s = Ant^a/m J n 2 (r)d 3 r, and the second one is 
due to the on-site dipole-dipole interaction Udd — 
l/(2w) J V(q)n 2 (q)d 3 q, being V(q) and n(q) the Fourier 



transform of the dipole potential and density, respec- 
tively [13[ . Due to the localisation of the wavefunctions 
at the bottom of the optical lattice wells, the long range 
part of the dipole-dipole interaction Ut is in a very good 
approximation given by the dipole-dipole interaction po- 
tential at distance i, U/= D 2 [l — 3 cos 2 (a jr)]/£ 3 , multi- 
plied by the densities and rij in the two sites. The 
quantity D is the dipole moment and ctg is the angle be- 
tween the orientation of the dipoles and £. 

The ratio between the total on-site interaction U = 
U s + Udd and the nearest neighbour dipolar interation 
Unn determines much of the physics of the system. It 
can be varied by tuning the on-site dipole-dipole inter- 
action Udd from negative to positive by changing the 
vertical confinement, or by changing the s-wave scatter- 
ing length via a Feshbach resonance, as recently demon- 
strated with Chromium atoms [14j ]. Alternatively one 
can consider to use heteronuclear molecules, or Rydberg 
atoms, which possess much larger dipole moments and 
will be hopefully soon available experimentally in optical 
lattices. 

The dipolar interaction potential decays as the inverse 
cubic power of the relative distance. In most theoreti- 
cal approaches the range is cut-off at certain neighbours. 
The precise choice of cut-off range and lattice size, de- 
termines the fractional character of the allowed ground 
state filling factors (e.g. for INN (or 2NN) only multiples 
of 1/2 (or 1/4) fillings are found). In the present work, 
we consider a range of interaction up to the 4th nearest 
neighbour and focus on the case of dipoles pointing per- 
pendicular to the plane of the lattice, where dipole-dipole 
interaction between atoms in the plane of the lattice be- 
comes isotropic, and in particular always repulsive. We 
consider the case of dipole-dipole interaction relatively 
weak (U/Unn — 20) and strong (U/Unn — 2) with re- 
spect to the on-site interaction [151 ]. 

We study the problem in the mean field regime, with 
an approach based on the Gutzwiller ansatz. This cor- 
responds to writing the wavefunction as a product over 
the different lattice sites (i) of single-site wavefunctions 
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In particular the time-depedence of the Gutzwiller coef- 
ficients fn allows to study the evolution of the state in 
real (t) and imaginary (r = —it) time 
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FIG. 1: (Color online) Phase diagram for weak and strong 
dipole-dipole interaction: U/Unn = 20 (a) and U/Unn = 2 
(b). The thick lines are the ground state lobes, found (for 
increasing chemicals potential) for filling factors equal to all 
multiples of 1/8. The thin lines of the same color are the 
metastable states at the same filling factor. The other lines 
are for filling factors equal to odd multiples of 1/8 [19|; some 
of the metastable configurations at filling factor 1/2 (I to III) 
and corresponding ground state (IV). Empty sites are light 
and site occupied with 1 atom are dark. 



The imaginary time evolution, which due to dissipa- 
tion is supposed to converge to the ground state of the 
system, in the presence of long-range interaction happens 
to converge often to different configurations, depending 
on the exact initial conditions. This is a clear sign of the 
existence of metastable states in the system. In the real 
time evolution, their stability manifests in typical small 
oscillations at frequency ojq around a local minimum of 
the energy. All the insulating metastable configurations 
present an insulating lobe in the J — fi phase space, as ex- 
plained below. They have a finite lifetime due to the tun- 
neling to different metastable states, which can be very 
long for small tunneling parameter J and large systems. 
Using a path integral approach in imaginary time 11611 . 
combined with a dynamical variational method (cf. [17|). 
we have estimated the tunneling time T to diverge for 
J — > and to scale like u>oT w exv[N s ft] exp[— N s hJ/ J] 
for J I J > 0.3, being N s the number of sites and J of the 
order of the tip of the lobe [18| . 

The most convenient method to determine the phase 
diagram of the metastable states is to use a mean-field 
approach perturbative in tp^. Performing the mean field 
decoupling of the Hamiltonian, the tunneling part at 
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FIG. 2: (Color online) Number of metastable states as a func- 
tion of /j, for weak and strong dipole-dipole interaction: (a) 
U/Unn = 20 and (b) U/Unn = 2. (c) Energy of the ground 
(thick line) and metastable states (thin lines) as function of u 
for strong dipole-dipole interaction (U/Unn = 2). The inset 
shows the energy levels at filling factor 1/2. 



1st order in the order parameter takes the form H t = 
— JJ2i('Pi a i + 'Pi a i)- Using the definition tpi = (a i e _/3ff ) 
in the limit (3 — ► oo, for a given classical configuration 
described by the density distribution m, one gets 
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where is the dipole-dipole interaction of one atom 
placed at site i with the rest of the lattice. A classical 
configuration rij is defined metastable , if there exists a 
region of the phase space J — fi (insulating lobe) , where 
Eq.(j4]) only allows the trivial solution (fi = 0, Vi. The in- 
sulating lobes exactly coincide with the stability regions 
found with the imaginary time approach. Studying the 
properties of Eq. , one can determine in a reliable way 
the insulating lobes for a huge number of configurations, 
which it would be impossible to access only by looking 
at the convergence of the imaginary time evolution. 

We investigate all possible configurations in 4 x 4 lat- 
tices with periodic boundary conditions, with all filling 
factors N a /N s (number of atoms/number of sites), rang- 
ing from 1/N S up to one, including the possibility of hav- 
ing double occupancy of the lattice sites. The quantities 
of interest that we extract from this analysis are: (i) the 
boundary of the lobes for the insulating configurations; 
(ii) for each value of the chemical potential, the num- 
ber of metastable insulating states present at very low 
tunneling; (iii) the energy of the ground state; (iv) the 
energy of all the insulating metastable states. 

Those results are summarised in Figs. Q]and[2] We ob- 



serve that for weak dipole-dipole interaction (similar to 
the hard core limit) the system presents an almost exact 
particle- hole duality (Figs 1 1121 a)). while for small on-site 
interaction, which allows double occupation of the lat- 
tice sites, many more configurations arise at filling factors 
larger than 1/2 (Figs lll2F b)). As shown in Fig[^c), there 
is usually a gap between the ground state and the lowest 
metastable state, which might allow to reach the ground 
state by ramping up the optical lattice under some adi- 
abaticity condition. However, this feature is strongly re- 
duced in the case of larger lattice sizes that we are going 
to discuss in the following. 

The number of metastable configurations and the vari- 
ety of their patterns increase very rapidly with the lattice 
size. Since for lattice sizes larger than 4 x 4 it is not pos- 
sible to track down systematically all existing configura- 
tions, we used a statistical approach where we run many 
times the imaginary time evolution for the same values 
of the parameters, each time changing the initial condi- 
tions. Exactly for the same reason why the metastable 
states exist, the convergence of such a procedure might 
be very very slow and is not always accurate. Hence, the 
stability of each of the obtained configurations is tested 
using the mean-field perturbative approach described in 
Eq.®, in order to confirm the existence of an insulating 
lobe. In general, we find that the configurations which 
differ from very regular ones by small defects are stable 
in a large region of the phase space, while the lobes cor- 
responding to configurations with many defects are very 
small. In FigO we show the insulating lobes for three 
configurations which differ from the checkerboard only 
by small defects. 

Very important issues are the initialisation and detec- 
tion of the atomic states in the lattice. One can use 
superlattices in order to prepare the atoms in configura- 
tions of preferential symmetry. This idea is pursued by 
several experimental groups [20] . We have checked that 
the presence of defects is strongly reduced when a local 
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FIG. 4: (Color online) Normalised spatial noise correlation 
patterns for configurations (I) to (III) in FiglTI [2^]. 



potential energy following desired patterns is added to 
the optical lattice. Note that the configurations obtained 
in such a way will also result stable once the superlatticc 
is removed, thanks to dipole-dipole interaction. 

The spatially modulated structures created in such a 
way can be detected via the measurement of the noise 
correlations of the expansion pictures [1, [2l|, I22T] : the 
ordered structures in the lattice give rise to different pat- 
terns in the spatial noise correlation function, equal to 
the modulus square of the Fourier transform of the den- 
sity distribution in the lattice. Such a measurement is in 
principle able to recognize the defects in the density dis- 
tribution, which could be exactly reconstructed starting 
from the patterns in the spatial noise correlation func- 
tion. The signal to noise ratio required for single defect 
recognition is beyond the present experimental possibili- 



ties. However, averaging over a finite number of different 
experimental runs producing the same spatial distribu- 
tion of atoms in the lattice, a good signal can be ob- 
tained. In FigfJJ we show the noise correlations for the 
metastable configurations at filling factor 1/2 shown in 
Figffl (I) to (III). 

The capability of initialise and read-out the state of 
the lattice makes those systems useful for applications 
as quantum memories. The controlled transfer of those 
systems from one configuration to another will be object 
of future studies. 

Alternatively to superlattices, structures obtained e.g. 
with atoms chips or microlenses arrays, where each lattice 
site can be addressed individually, could be used to pre- 
pare the desired configuration and manipulate it. Finally, 
it will be worth investigating the possibility of creating an 
atom-light interface to initialize or read-out the atomic 
state of the system by coupling it with the polarisation 
degrees of freedom of light 24], for instance exploiting 
the spinor character of Chromium dipolar atoms [l2l ]. 
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